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Abstract 

A quantum mechanical picture, relating accelerated geodesic deviation to creation of massive particles via 
quantum tunneling in curved background spacetimes, is presented. The effect is analogous to pair production 
by an electric field and leads naturally to production of massive particles in de Sitter and superluminal FRW 
spacetimes. The probability of particle production in de Sitter space per unit volume and time is computed 
in a leading semiclassical approximation and shown to coincide with the previously obtained expression. 



1. When embedded into an external background, the "vacuum" of quantum field theory can lose stability 
because of particle creation. The difference from classical or quantum radiation is that the particles are produced 
through quantum tunneling. 

The standard paradigm is the "Klein paradox" [T], i.e. creation of electron-positron pairs by an electric 
field E. Let us for simplicity consider a constant electric field along the z— axis in a capacitor of width L. The 
probability to create a single particle-antiparticle pair (the Schwinger process [2]) in its center of mass frame is 

i 2 

given by the tunneling exponential exp(i J p(z) dz)\ , where the momentum p(z) is the imaginary solution of 



the dispersion relation 

(S+ -eAo(z)) 2 = 4p 2 + 4m 2 , (f) 

with the total energy of the pair £ + = —2m. In the above equation z = z\ — Z2 is the distance of two particles, 
p = (pi — P2)/2 is its conjugate momentum, while Ao(z) = J z Edz. 

The characteristic parameter is the "creation length" Z, defined by the condition that the energy offered by 
the electric field compensates for the rest energy on the particles produced, i.e. eAg(l) = e J Q Edz = 2m, which 
gives I = 2m/ eE for a constant field E. The pair creation occurs only if L > I and the probability of creation 
per unit volume and unit time interval is the well known exponential [5] 

6{L - I) cxp (- yjArn 2 - e 2 E 2 z 2 dzj = 9(L - \)e-^l eE . (2) 

This integral runs over the distance z between the particles in the pair and describes the tunneling from the 
point z = at the "moment" of creation to z = I. Note that the applicability of the quasiclassical approximation 
requires the exponent in ^) to be large and negative, so that the exponential is small. 



2. The situation becomes considerably more interesting when the electric field is replaced by a gravitational 
one [3]. For instance, contrary to the case of a constant electric field, which can pull apart the two virtual 
particles of the pair, the equivalence principle forbids pair creation in a constant gravitational field. 
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1 The energies E\, £2 and momenta pi,P2 along the z— axis of the particles in the pair are related by 



■ eEz-)±- (£- 
' 2 v 



■ eEz+) 
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(P+ ±P-) 



with z± = z\ ± 22, £± = £1 ± £2, p± 



which in the center of mass frame p+ 



- pi ± p2 ■ Adding these two equations one obtains 
{£+ - eEz-) 2 + (£- - eEz+) 2 

- 0, 2+ = and £— = 0, reduces to JTJ 



p\ +pL + 4m 2 . 



f 



However, it is natural to expect that an analogous "pulling-apart force" will occur also in gravity, in the 
presence of curvature, as long as the latter has the right properties, which make close geodesies diverge. 
The equation of geodesic deviation in a given background specetime is [3] 

D\rf = R\ pa u v u p rf , (3) 

where R^ vpc , is the Riemann curvature tensor, w M = dx^/dr is the four velocity along the "central geodesic" 
of the one-parameter congruence x^{t, A), parametrized by A, and rf = dx^/dX, is essentially the coordinate 
distance between two close probe particles moving along the corresponding geodesies. D T denotes the covariant 
differentiation with respect to the afhne parameter r along the central trajectory. 

Consider the case of a maximally symmetric background spacetime with constant curvature scalar R. De 
Sitter and AdS are two obvious examples of interest. The Riemann tensor of such a space is of the form 

R 

R^vpa — — Y2 (.9n<?9vp — 9^p9v<y) (4) 

Inserting Q into ([3]), one obtains 

0V = ~ fa" - Mu") . (5) 

In the study of the particle creation process, in particular, one can choose if such that {urj) = 0. That is, in 
the rest frame of one of the particles with u = (1,0, 0,0), one observes the motion of the second one at a proper 
space-like distance from the first, i.e. with 77 = 0. Then, equation ([5]) simplifies to 

DW = -^V\ (6) 

where Latin indices denote the spatial coordinates. 

In eqs.©-© can not be chosen arbitrarily: it is a tangent vector to a family of geodesies and is related 
by the zero-curvature condition to u^. Physically significant is, however, not rf" , but its length z = \\r)\\, which 
measures the distance between the geodesies. In any spatially isotropic spacetime one may choose rf to have 
only one non-zero component. In this case, using (|6]) one obtains for the equation of motion of : 



d 2 z R 

This equation, in turn, can be converted into the dispersion relation 



(7) 



where p is the momentum conjugate to z. It is amusing to notice, that this equation describes a non-relativistic 
inverted harmonic oscillator and the tunneling probability is again given by the integral between the turning 
points and is proportional to (after putting £ — —2m) 



cxp^-2^ = 0{-R)exp ^- 2 ^ 



Am+ _| Z 2 dz \ = exp \- J™ 12 ) ■ (9) 



2 Indeed, 



D 2 T z 2 = 2g^n li D 2 .ri v + 2 gt!lv D T rf D T rf = g^rfjf + Ig^D^rf D T rf = z 2 + 2g^D T TfD T rf . 
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On the other hand, 



Therefore, 



„9 9 & 2 Z 



d 2 z R 

K = Z+X, 

dr 2 12 



where 



X = -^-^D^^Drrf - z 2 g^D T ^ D T rt") . 



1 

For any spacelike vector 7/ M with only one non-zero component X = 0. 

A simple example to illustrate these points is S 2 . For meridians <f> = const on a two-dimensional sphere with the metric 
dB 2 + sin 2 8d<j> 2 the role of afhne parameter r is played by 9, and the two tangent vectors are u = (1, 0) and n = (0, 1), while 
z = WvW = sin 6 and D^r] = —n while d 2 z = —z. 
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Again, the applicability of the quasiclassical approximation used here, requires that the exponent in this formula 
is large, i.e. m 2 3> —R/12 = H 2 . 

This formula implies that particles are produced in de Sitter space when the Ricci scalar is negative, 
R = —12H 2 < (H is the Hubble constant) with probability per unit volume and time interval given by 
exp {—2irm/H) [5], while they are not produced in AdS space, where the Ricci scalar has the opposite sign and 
the relevant geodesies decelerate rather than accelerate □ . In the case of massless gravitons in de Sitter, the 
role of m 2 is played by k 2 /a 2 (t), for a given mode with comoving momentum k. The production of modes with 
k 2 /a 2 (t) < H 2 is expected not to be exponentially suppressed [5]. 



3. Similarly, one can deal with more complicated non-maximally symmetric spacetimes. In order to de- 
termine if particles are produced in these cases, it is not sufficient to look at the sign of the curvature scalar; 
instead, one has to proceed with a more refined analysis of the geodesic deviation equation. For instance, 
consider the Friedmann-Robertson- Walker (FRW) metric 

ds 2 = n 2 (x°) ((dx ) 2 - (dx) 2 ) = dt 2 - a 2 (t)(dx) 2 (10) 

and assume for simplicity that the 3-dimensional space is flat. Then, the curvature scalar is equal to 

il" 



where the prime denotes differentiation with respect to x , and the geodesic deviation equation reads 

= (%■- 2^ ) ((«°)V - (A°K) + % (v l - (mW) . (12) 



In the case of de Sitter, n(x°) = l/(Hx°) = exp(Ht). Thus, the first term in this equation vanishes and one 
returns to formula ([5]). In the generic FRW case, on the other hand, in the comoving coordinates of the central 
geodesic i.e. for i// 1 = (1,0,0,0), equation (fT2"j) takes the form 

-(£-£)''-!<'■ (13) 

where the dots denote differentiation with respect to t. Again, this equation can be rewritten as an equation 
for the distance z between two neighboring geodesies, with the covariant derivatives replaced by the ordinary 
ones with respect to t, namely 

d 2 z a , , 

W = a Z ' (W) 
This equation also describes an oscillator but with time-dependent frequency, and, depending on the explicit 
form of a(i), there can be relatively accelerating geodesies, in which case particle creation will take place. 

The lesson is that the particle production is determined in general not simply by the sign of the curvature 
scalar, but by the sign of d (given that a > 0). Pairs are produced only if a > 0. Note that a solution of 
(fT4|) is just z(t) = a(t). For a power-law behavior of the scaling factor a(t) ~ t s , the condition d > becomes 
s(s — 1) > 0. This inequality is satisfied for s > 1 or s < 0. The latter behaviour, however, describes a 
contracting universe in which, of course, there is no chance of particle creation. On the contrary, in the case 
of s > 1 describing a superluminal accelerated expansion, pair production takes place. In order to calculate 
the production rate in these FRW metrics, one has to deal with an oscillator with varying frequency and, 
correspondingly, a time dependent barrier. This requires detailed analysis and the answer will depend on the 
explicit form of a(t) [7]. 



4. A few simple thoughts related to the fate of de Sitter evolution as a result of particle production are in 
order. Since emerging particles gravitate, they act against the "anti-gravity" of de Sitter space, which "repels" 

3 The emergence of the oscillator potential JHJ in the problem related to de Sitter space should not come as a surprise: in the 
static coordinate system the de Sitter metric has the form 

ds 2 = (1 - H 2 r 2 )dt 2 ^—r^dr 2 ~ r 2 dQ,l, 

1 — H 2 r 2 

which in the Newtonian approximation corresponds to the inverted oscillator potential V = —H 2 r 2 /2. 
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particles. It is an old hypothesis [5], that maybe the back reaction of created particles could drastically affect 
the fate of the de Sitter evolution, perhaps even stopping its accelerated expansion. For this to happen, however, 
the density of created particles should be big enough and maybe "explosive" . This is not easy to achieve due 
to the dilution of the produced particles by the space expansion. 

The particle density n — N/V in the absence of self-interactions is governed by the well-known equation 

V 

h = a — n— , (15) 

where the first term on the r.h.s. describes the creation of particles, a ~ exp (— 2-Km/H), while the second 
term describes their dilution. In the exponentially expanding flat de Sitter one obtains V ~ exp(3_ffi), so that 
V/V = 3H = constant and (fT5|) implies that 

n=^+/3e" 3 ™, (16) 

which does not grow with time and the first term is exponentially small so that the creation of particles cannot 
stop the accelerated expansion. 

This well-known result is due to the exponential expansion of de Sitter space. One way one might expect 
to change the situation is by introducing self-interaction among the produced particles (see [9j [10] for further 
discussion) . The self-interaction can stimulate production of particles with a chain reaction, represented in (|15p 
by an extra term of the form 771 with 7 > on the right hand side. Hence, if 7 > 3H , this term dominates 
over dilution and particle production can eventually cancel the de Sitter exponential expansion, as a result of 
the exponential increase of the matter density 0: 

+ (3e^- 3H » . (18) 
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No quantum correlations, fusion of created particles, or the effect of Bose statistics was taken into account 
above. The latter, in particular may lead to higher powers of n on the r.h.s. of the density evolution equation 
and has been discussed in [3]. 

5. Any account of self-interaction which leads to fast increase of particle density implies that the Einstein 
equations should be modified. In particular, this means that the acceleration of the Universe expansion can 
change to deceleration. Indeed, 

-~-(3p + e) (19) 
a 

where e is the energy density and p is the pressure. Since in the de Sitter space e = H 2 — —p, the acceleration 
is positive. However, if enough particles are produced with positive both e and p (the relation between e and p 
depends on the model), one can change the sign of the r.h.s. of p^|) . 
However, the curvature 

fl~-(e-3p) = -T" /4 (20) 

only increases!! in its absolute value when any real particles are created, with e — 3p > 0. Thus, while creation 
of particles can indeed act against and, perhaps, even stop the accelerated expansion, looking at the curvature 
and assuming that it is the negative curvature which controls particle production one could, at first sight, think 
that there is violation of the Le Chatelier principle, namely, that negative curvature creates particles, and at the 
same time it becomes more negative as a result of this. However, as we saw in Section 3, it is the acceleration a 
and not the curvature that controls the particle production. Since, according to (|19[) the acceleration decreases 
as a result of particle creation, there is no contradiction with the Le Chatelier principle. 

6. An elementary quantum mechanical (as opposed to quantum field theoretic) discussion of particle pro- 
duction in gravitational backgrounds was presented, appealing as much as possible to physical intuition and 



4 If the spatial section of the space-time is closed, i.e. V ~ cosh 3 (_fft), one obtains instead 



S costr 3 (Ht) 



3ff-7 3H + -y H--y + 7 



■ fi Jm> (17) 



which does not change the conclusion. 

5 I n |H| . the authors obtained a decrease of the cosmological constant as a result of negative corrections to e due to the interaction 
between gravitons. 
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carefully avoiding the formal and usually technically very involved argumentation. It relates the phenomenon of 
pair production to the accelerated deviation of nearby geodesies in the given spacetime background. In simple 
cases of maximally symmetric backgrounds, it leads to the known exponential barrier penetration formula for 
the probability of such particle production. The approach is general, valid in locally curved spacetimes of any 
dimensionality. Our arguments are not expected to be new to the experts, but surprisingly, they are not, to the 
best of our knowledge, explicitly presented in the existing literature. 
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